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ABSTRACT 

Using 2-loop renormalisation group calculations, we study a system of N two- 
dimensional Potts models with random bonds coupled together by their local energy 
density. This model can be seen as a generalization of the random Ashkin- Teller 
model. We found that, depending on the sign of the coupling term, the universality 
class of the system in the presence of randomness is different. Under particular con- 
sideration, this model presents an example of a first order phase transition rounded 
by randomness. 
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The effect of quenched randomness in two-dimensional systems having a continuous phase 
transition in its absence has been studied for a large variety of cases. In particular, using a 
renormalisation group (R.G.) approach, the effect of uncorrelated impurities has been studied 
for the two-dimensional Ising model Baxter model and for the Potts model |4j. For 
most of these cases, the Harris criterion || provide us a good method to see if disorder will 
change or not the universality class. For first order phase transition systems, the effect of 
randomness can be stronger. It was argued in Q and established in a rigorous way in 0] that 
for two-dimensional systems the transition becomes continuous in the presence of impurities. 
Recently, it was shown by Cardy § that such phenomenon effectively happens for a large class 
of systems presenting a fluctuation driven first order transition in the pure case. In view of 
the results obtained there, it was suggested that the Ising-like transition could be a general 
feature of two-dimensional systems with randomness. Of particular interest is the random N 
color Ashkin- Teller model considered in this work. This model can be easily generalized to 
N coupled g-state Potts model (with 2 < q < 4) with random bonds. The aim of this paper 
is to compare the behavior of this model in the pure case and in the presence of quenched 
randomness for generic N. Our result is that when impurities are added, the system will 
present a continuous transition but with different universality classes depending on the sign of 
the coupling between the models. In any case, the transition is not Ising-like. 

Our model consists of N two-dimensional g-state Potts models near the critical point coupled 
together by their energy operators. This model can be considered as a generalization of the N 
color Ashkin- Teller model with 2 < q < 4. The Hamiltonian of the system has the following 
form: 

N . N . N 

H = ^H 0)i + / d 2 x m^£j(x) - g d 2 x^£ i £ j {x) . (1) 

i=l ' i=l ' iy^j 

i?o,i are the Hamiltonians of the N Potts models at the critical temperature, m is the reduced 
temperature, e corresponds to the energy operator of the pure models and the last term in (||) is 
the coupling between the different models. This model was well studied in the case q = 2 which 
is the N color Ashkin- Teller model. For N = 2 it turns to be integrable ||. For N > 2 the large 
scale behavior depends on the sign of g: For g < it presents an Ising type second order phase 



transition while for g > a fluctuation driven first order phase transition is present [10, 11]. 



Recently, it was shown [12] that for N = 2, q > 2 the model is still integrable presenting now 



a mass generation indicating again a first order phase transition. The situation is not so clear 
for N > 2. In the lack of an exact solution we can give a R.G. approach by studying the large 
scale behavior of the system for different values of g. This was done in || |j] for the case N — ► 
in the context of the Potts model with random bonds. A non trivial fixed point was found with 
critical exponents for the specific heat and magnetic susceptibility which differ from the ones 
of the pure model. 

For generic N, in the limit m — > 0, the 2-loop R.G. equations of (Q) can be obtained from 

I: 

P(g) = = eg(r) + 4n(N - 2)g 2 (r) - 16vr 2 (iV - 2) 5 3 (r) + 0(g\r)) 
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Here, Z e and Z a are respectively the renormalisation constants of the energy and spin operators 
and e = 2 — 2A e where A e is the physical dimension of the energy operator of the pure model 
(e differs from that of Q by a factor of —3, so e = for q = 2 and e = | for (7 = 3). For the case 
iV > 2, we can see that if g > the coupling constant flow far from our perturbative region. 
Even if a definite proof is not given, a comparison with the case q = 2 seems to tell us that a 
mass gap is dynamically generated indicating a first order phase transition. The situation is 
completely different for g < 0. In this case, we can see that there is a non trivial infrared fixed 
point at: 



c 2 



^ = -4^V^) + 4 vr(iV-2) 2 +Q(£) (3) 

The critical exponents associated with the energy and spin operators in this new infrared fixed 
point are given by: 

a; = a„ - 7c ( So ) 

= A, + g41 £ - MzQj + 0l f) (4) 



and 



(N - 2) 2(N - 2f 
A' a = A a - la (g c ) 

( (N-i) \ r 2 (-|)r 2 (i; 



1 



+ 0(e 4 ) (5) 



a V.8(iv-2)V r 2 (-|)r 2 (- 6 , 

The results obtained in [||, [|] for the case of the quenched random case can be reproduced by 
just putting N = in the formula above and keeping in mind that these results then make sense 
only for g > 0, which in this case g corresponds to the variance of the probability distribution 
for the random bonds. 

Let's consider now adding quenched randomness which couples to the energy operator (e.g. 
the local energy density). This can be done by introducing in (|l]) a position dependent random 
mass term m — > m{x) where m{x) = and m(x)m(y) = A5(x — y), with A > 0. Using 
the replica method, we can introduce n copies of the system and average over a Gaussian 
distribution for m, this will give us a Hamiltonian of nN Potts models coupled by their energy 
operators. The replicated Hamiltonian is: 

H = Y.HS, t -gjd 2 xJ2^(x)-^Jd 2 x £ ete)(x) (6) 

where indices i, j runs from 1 to N and a, b from 1 to n and < • • • > means that terms with same 
replica indices are omitted. As explained in H, this is because these terms produce irrelevant 
operators or trivial contributions which could just give a shift in the critical temperature. This 
model was studied at 1-loop level in || for the case q = 2 (the random Ashkin- Teller model). 
It was shown that when randomness is present, the R.G. trajectories for the coupling constants 
run away in a first stage from the U.V. fixed point but finally curl around and approach the 
fixed point corresponding to iV decoupled Ising models. 

The procedure to obtain the 2-loops R.G. equations for @ is the same as in 0]. After some 
combinatory, we find, taking directly the limit n — > 0: 

g = eg + (N - 2)g 2 - 2gA - (N - 2)g 3 - (2N - h)g 2 A + 4gA 2 

A = eA - 2A 2 + 2(JV - l)Ag + 2A 3 — (N — l)g 2 A - 2(N - l)A 2 g (7) 
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where for simplicity we have made the change g — > Airg, A — > 4ttA. The first step in studying 
the flow given by (Q) is to look for the fixed points g = A = 0. Let's first consider the case 
N > 2. Then, we find four fixed points given by: 

g = ; A = (8) 
+ ,*\,o +0(e 3 ) ; A = (9) 



(JV-2) (iV-2) 



e 2 



g = ■ A = - + - + 0{e i ) (10) 

The others fixed points present in ([?]) are of order of constant and so beyond our perturbative 
region. The non trivial fixed points (g) and ([1C|) are the ones present in the case of iV non 
random coupled models and one random model respectively. There is however a new fixed 
point ([ll]) that can not be seen at the 1-loop level. As usual, we can study the stability of each 
of these fixed points by re-expressing (0) around the solutions above g = g c + 5g, A = A c + 5 A 
and keeping only the smallest order in e. This will give us a linear system: 

5a)- A \5A 

The information about the stability of each of these fixed points can be obtained by calculating 
the eigenvalues of A in each of these cases. In this way it is easy to see that (||) is unstable (for 
e > 0), (0) and (|ll|) are both stable while fli~0| ) is stable only when g = 0. This result tell us 
that the random fixed point found in j^, |j can not be reached if a small coupling between the 
different Potts model is added. To study the coupling constant flow , let's first suppose that 
the initial conditions are A(0) , g(Q) > 0. It is useful in a first stage to keep just the first and 
second order terms in (0). In this case, it can be shown that for generic e the solution is given 
by: 

A(0 = erA o (r) ; g(l) = erg Q {r) (12) 

where r = e d /e, I is the R.G. scale parameter and Ao, go are the solutions of (0) with e = 0. 
A closed form for Ao and go has been given in ||, it is easy to see that for r — > oo we have 
go (t) << Ao(r) and Aq(t) ~ So, the trajectories of the coupling constants for generic e will 
flow toward the region: g(l) « A(Z) and A(/) ~ |. This is precisely the region where third 
order terms in ([/]) become important indicating that trajectories will asymptotically reach the 



point (11). We can now calculate the critical exponents associated with the energy and spin 
operator. Since the difference between the points ( JlOj ) and (11) is of order e 2 we can not see at 
this order in perturbation any difference in the spin critical exponent compared to that found 
in [U for one Potts model with random bonds. For the energy operator, the correction to the 
critical exponent up to second order is given by: 

A' £ = A £ + A C -(N- l)g c - ^ + ' • • 

= A £ + | + e j + 0(6 3 ) (13) 

giving surprisingly for generic ./V the same result as for N = 1. Repeating the same arguments, 
we can show that for initial conditions A(0) > , g(0) < the R.G. trajectories will now flow 
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toward the point (||) which is the nonrandom fixed point. The critical exponents will be given 
by (|^) and (H|) indicating that for this sign of g randomness will not change the large scale 
behavior of the system. 

Let's finally turn to the particular case N = 2. As was said before, the non random case 



is exactly integrable [12] giving a mass for any non zero value of g. It is interesting to see if 
randomness will increase the order of the transition or will keep a non zero mass in the theory. 
The R.G. equations are simply obtained by rewriting (0) for N = 2: 

g = eg-2gA + g 2 A + 4,gA 2 

A = eA - 2A 2 + 2Ag + 2A 3 - g 2 A - 2A 2 g (14) 



The fixed points (g|) and (|l0[) are still presents with the same type of stability and ( Jll[) becomes 
now: 

; A = i±^ (15) 

But now, as expected, (§) has disappeared. In the region g > the arguments given for generic 



iV are still valid and trajectories will flow toward the point (15) giving at this order the same 
spin and energy exponents as for generic N. The difference appear for initial conditions with 
g < 0. In this case, the R.G. trajectories will flow far from our perturbative region. With 
the results above we can not say if trajectories will flow toward a massive theory or a non 
perturbative fixed point. A naive analysis of (0) in this case indicates that trajectories will 
reach asymptotically the integrable line A = which gives a massive theory. This would be 
in contradiction with Imry and Wortis arguments || and statements given in 0. A similar 
situation happens in the case q = 2, N = 2 (the random Baxter model) studied in [||. This 
example of a system with quenched randomness violate the Harris criterion which is in some 
sense the second order transition version of Imry- Wortis argument. 

In this paper we have analyzed the behavior in presence of quenched randomness of a model 
which present a very rich structure. Depending on the sign of the coupling g, the nature of the 
transition is different in the pure case. We found that in the random case, for N > 2, whatever 
the sign of g is, the transition is continuous and different from that of the random Ising model, 
but the universality classes are different for g > and g < 0. Moreover, in the case g < 
the disorder doesn't change the universality class of the model. This apparent contradiction of 
Harris criterion is due to the fact that randomness, which is relevant at the tree level, is driven 
to be irrelevant by the coupling term. The behavior of the model in the region g > is also 
surprising: the fixed point of N > 2 coupled Potts models in the presence of impurities is not 
the same as found in [||, |I| in the case N = 1. The spin and energy exponents remains however, 
at this order in perturbation, the same for any value of N. 

The particular case TV = 2 is more complicated; for g > the behavior of the R.G. flow is 
similar to the one of the case N > 2. Since the pure case is integrable and is a massive theory, 
we have a system which certainly gives a first order phase transition in the pure case. The 
presence of randomness makes however the transition continuous. This give a new example of 
a first order transition rounded by randomness. It is not clear however what finally happens 
for g < and only a non-perturbative treatment would give a definitive answer. 

I would like to thank VI. S. Dotsenko and M. Picco for many helpful suggestions and 
conversations. 
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